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Abstract 

We numerically demonstrate that "mode-entangled states" based on the transverse modes of 
classical optical fields in multimode waveguides violate Bell's inequality. Numerically simulating the 
correlation measurement scheme of Bell's inequality, we obtain the normalized correlation functions 
of the intensity fluctuations for the two entangled classical fields. By using the correlation functions, 
the maximum violations of Bell's inequality are obtained. This implies that the two classical fields 
in the mode-entangled states, although spatially separated, present a nonlocal correlation. 
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Quantum entanglement is a fundamental concept and one of the most interesting prop- 
erties of quantum mechanics. The importance of quantum entanglement in quantum infor- 
mation processing is by now widely appreciated [l|, Quantum entanglement is first 
introduced by Einstein, Podolsky, and Rosen (EPR) as most noticeable the EPR paradox 



4] 



which is at the origin of quantum nonlocality. Bell proposed a remarkable inequality 
imposed by a local hidden variable theory, which enables an experimental test on the quan- 
tum nonlocality [5j|. The violation of Bell's inequality indicates the presence of quantum 
entanglement. Generally, the nonlocal properties, which appear in the correlation measure- 
ment of quantum entangled states for two spatially separated particles, are considered as the 
inherent features of quantum theory with no classical analog. Numerous theoretical studies 
and experimental demonstrations have been carried out to understand the nonlocal proper- 
ties of quantum states. It is noticeable that the violation of Bell's inequality is possible in 
nonlinear optical processes such as multimode parametric amplifiers and four-wave mixing, 
in which the output modes are in some way correlated jo) . In order to explore the intrinsic of 
quantum entanglement, some researchers yet use classical optical fields to simulate quantum 
entanglement numerically and experimentally 0, El- 

The similarities between the Helmholtz equation and the Schrodinger equation have at- 
tracted some researches on the analogies between the transverse modes in multimode waveg- 

1(1 11 1 . Besides the uncertainty relation jl^j], Wigner 



uides and the quantum Fock states [9 , 



phase-space distributions of the optical multimode fields exhibiting negative regions are sim- 
ilar to quantum Fock states, which had been verified experimentally ^| . Farther, some 
researchers have argued convincingly that the classical Maxwell field plays the role of the 
quantum wave function for a single photon Q|. This might be the physical basis of the 
similarities. However, the similarities have been limited principally to the measurement of 
first-order coherence, i.e., single-particle states. Classical wave analogs of the high-order non- 
ocal coherence (quantum entanglement), i.e., multiparticle states, have been seldom studied 
3]. The researches on the high-order local coherence can date back to the theoretical and ex- 
perimental researches on interferences of two independent laser beams in Refs. ^,12]. Due 
to the random phase difference between the two beams, no stable but transient interference 
fringes can be observed. These local effects can be well explained by quantum theory and 
classical electromagnetics. Recently, "mode-entangled states" based on the transverse modes 
of classical optical fields propagating in multimode waveguides are proposed as the classical 
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simulation of quantum entangled states [3]. The states can be regarded as the nonlocal 
generalization of the high-order local coherence. It is interesting that the mode-entangled 
states can also show the nonlocal correlations such as the violation of Bell's inequality and 
the nonlocal properties of optical pulses' group delays. However, the nonlocal correlations 
are generally considered as the inherent features of quantum entanglement. 

In this letter, we numerically demonstrate that the mode-entangled states violate Bell's 
inequality by completely using classical electromagnetics. Numerically simulating the cor- 
relation measurement scheme proposed in Ref. |l8j ]. we obtain the normalized correlation 
functions of the intensity fluctuations for two entangled classical fields. Then the correlation 
functions are substituted into Bell's inequality and the maximum violations of Bell's inequal- 
ity are obtained. Despite all that, the mode-entangled states and the quantum entangled 
states, from a physical viewpoint, are different because the measurements of classical fields 
and quantum states are different. The research on this simulation may be important, for it 
not only helps to understand the nonlocal properties of quantum entanglement from a new 



viewpoint, but also arouses interest in a full 



the transverse modes of classical fields 
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optical quantum computation scheme based on 



201 ] . And besides, it may be attractive to other 



research fields. Such as researches on matter waves of Bose-Einstein condensates, similar to 
the scheme, a quantum computation scheme based on the transverse modes of the guided 



matter waves 211 is also interesting. 

n 

In the correlation measurement scheme of Bell's inequality proposed in Ref. 18], two 
independent monochromatic classical fields are prepared to input a CNOT gate as the control 
and the target fields, respectively. When the control field is given at mode superposition 
states (|TE ) ± |TEi)) /v2 and the target field is given at the mode |TE ) or |TEi), the 
output states of the CNOT gate are so called "mode-entangled states" , 

|d>±) = -L (|TE ) C |TE ), ± |TEx) c |TE0 t ) , (1) 



\9f) = 4 (l TE o) c |TEi) t ± |TE X ) C |TEo) t ) , 



V2 
1 

where subscripts c and t represent the control and the target fields, respectively. The states 
in each waveguide are mode superpositions, which are obviously different from product 
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states, 

|* 2 ) = |* c ) ® \%) (2) 
= ^(|TE > c ±|TE 1 ) c )(|TE > t ±|TE 1 > t ), 

where |^ c ) and \^ t ) represent the states of control and target fields, respectively. The differ- 
ence can be obtained not by measuring the single field, but by the correlation measurement 
of the control and the target fields. The correlation measurement can present that the two 
entangled fields are inseparable to some extent, which is a nonlocal correlation. 

To perform the correlation measurement of Bell's inequality, the control and the target 
fields output from the CNOT gate are sent to spatially separated (independent) mode an- 
alyzers (MAs), each of which contains a Y splitter and a variable phase modulator 6\ (9 2 ). 
Based on the correlation analysis, we get the correlation function, 

s(e 1 ,e 2 ) = (A c (e 1 )B t (e 2 )), (3) 

where A c (9i) and B t (9 2 ) are the intensity difference operators of the MAs' outputs for the 
control and target fields, respectively, 

A c (9 1 )=i^9 1 )-i;(9 1 ) (4) 
= e 2 ^ |TE ) C (TEx| c + e^ 1 |TEx) c (TE | C , 

Bt (9 2 ) = /+ (9 2 ) - I- (9 2 ) 

= e 2ie * |TE„) t (TE^ + e- 2i02 \TE 1 ) t (TE | t , 

with if (9i) and if (9 2 ) being the MA's operations on the control and the target fields, 
respectively. Apparently, there is no correlation between the independent measurements of 
the two fields in the product states, hence, their correlation function can be written as the 
product of the expected values, 

S* 2 9 2 ) = (* 2 | A c (9 1 ) B t (9 2 ) |* 2 ) (5) 
= (*c|i c (0i)|*c)-(*t|A (6 2 )\%). 

However, the correlation function of the mode-entangled states can not be written in such 
form because the results of the independent measurements of the two entangled fields are 
correlated. Here Bell's inequality is given as the criterion to distinguish the mode-entangled 
state and the product state. 
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Consider an ensemble (large collection) of the independent and identical systems that 
are labeled by A that satisfies normalization condition / p (A) dX = 1, where p (A) is a 

J A 

distribution function and A is spanned by A. In this model, the correlation function can be 
rewritten as a normalized form 

S(9 1 ,9 2 ) = (A c (9 1 )B t (9 2 )) (6) 
A' c (6*i, A) B' t (9 2 , A) p (A) dX, 



with 



A' c (0i,\) = m , (7) 

[J A ^(^,A)p(A)dA] 1/2 

B' t (9 2 , A) - Bt{ ° 2A) 



[J A BUe2,x)pWdx] 1/2, 

where A c (9\, A) and B t (9 2 , A) are the correlation measurement results corresponding to each 
A. Using the property of the product state's correlation function shown in Eq. Schwarz 
inequality \S (9i,9 2 ) \ < 1 and Eq. ©, we obtain 

\S(9 U 9 2 ) -S(9 1 ,9' 2 )\ <1± / A> c (9[,X)B>(9' 2 ,X)p(X)dX (8) 

J A 

+ 1± / A' c (9[,X)B' t (9 2 ,X)p(X)dX 

J A 

= 2±[S (9' 1 ,9' 2 ) + S(9' 1 ,9 2 )}, 

then Bell's inequality (CHSH inequality j^]) is given explicitly as 

\B\ = \S(9 1 ,9 2 )-S (9 1 , 9' 2 ) + S (6[, 9' 2 ) + S (6[, 9 2 )\ < 2, (9) 

which implies that the violation of Bell's inequality never occurs for the product states. On 
the contrary, if the violation occurs, there is a nonlocal correlation between the two fields 

n 

of the mode-entangled states. In Ref. p, Bell's inequalities are generalized for the case of 
fields or arbitrary intensity and the violation of Bell's inequalities is shown to be possible in 
a regime showing strong violation of Schwarz inequality. But the Schwarz inequality, related 
to the phenomena of photon antibunching and squeezing, is different from \S (9%, 9 2 )\ < 1 
satisfied by the correlation function. 

Let us now numerically simulate the scheme shown in Fig. 1 to obtain the correlation 
functions of the control and the target fields. For such purpose, we have to simulate the two 



independent optical fields propagating in waveguides. In Refs. jlall2|) the two independent 
light beams are replaced by two monochromatic optical fields with random phases <p\ and (fi2 
uniformly distributed in [0, 2n]. Therefore, the two independent optical fields propagating in 
the waveguides can be regarded as an ensemble labeled by the phase difference A = tpi — (f2, 
which also is a random variable uniformly distributed in [0, 2ir]. Thus, a simulation of each 
A is equivalent to an independent experiment of the ensemble, then the integrals of the 
simulations describe the behaviors of two independent optical fields. Although it requires 
infinite numbers of the simulations in theory, we can assign A with discrete values in practice 
to realize the ergodicity of A by finite numbers of the simulations. 

The finite differential beam propagation method (FD-BPM) employed in our numerical 
simulation is widely used in the analysis of optical waveguides and well proved effective and 
reliable by a lot of researchers |23j . Utilizing FD-BPM, the evolution of waveguide modes 
are clearly demonstrated using intensity distributions. And the waveguide output intensities 
can be obtained by the integral of the intensity distributions. In the simulation of the scheme 
in Fig. 1, the two MAs are separated far enough to avoid disturbing each other, and the 
numerical simulation result is shown in Fig. 2. For a certain A, we obtain the relation 
between the MAs' output intensities and 8\ (6 I 2 ) by changing the refractive indexes of the 
phase modulators. We find that the output intensities of MAs also change with the phase 
difference A. The reason of the fact is that TE modes of the control field and the target field 
in one of the Mach-Zehnder interferometer's arms, due to the phase difference A, interfere 
with each other when they are passing through the directional couplers of the CNOT gate. 
The interference causes the intensity in Kerr-like medium varying with A, therefore the phase 
shift via the Kerr effect also varies. Finally, the MAs' output intensities change with A due 
to the phase shift. The intensity differences A c (9i,X) and B t (02,X) versus 6\ and 62 are 
shown in Fig. 3. While A varies in [0, 2ir], A C (6\,X) and B t (6 2 ,X) fluctuate around their 
mean values, and the ranges of the fluctuations change with 6\ and 62, respectively. 

First, we obtain the random sequences A c (6\, A) and B t (62, A) subjected to a sufficiently 
long random sequence of A, then substitute them into Eq. (JHJ), and obtain the correlation 
function S (61, 6 2 ) via summation of A c (6±, A) and B t (82, A) instead of the integral of A. The 
correlation functions S (61,62) of and are shown in Fig. 4. After substituting 

S (61,62) into Eq. we obtain the maximum violation of Bell's inequality, as shown in 
Table 1, where the maximum values of \B\ are the average results of many A's sequences. 



6 



Considerable attention should be paid to that S {61,62) shown in Fig. 4 are quite different 
from the correlation functions cos (26 1 ± 262) or sin (26 \ ± 262) of the quantum entangled 
states, therefore, 6\, 6[, 62 and 6' 2 of the maximum violation of Bell's inequality in Table 1 
may not be 7r/8, — vt/8, 0, and 7r/4. During the simulation, we find that the granularity of 
A hardly influences the simulation results, so that the A is assigned with 64 discrete values 
in [0,27r]. Farther, we obtain the relation between the maximum of \B\ and the refractive 
index of the phase modulator 6 3 , as shown in Fig. 5. 

In this letter, we have numerically demonstrated that the mode-entangled states gener- 
ated by two independent classical fields propagating through the CNOT gate are different 
from the product states. The difference is presented in the correlation measurement of the 
two fields. For the mode-entangled states, Bell's inequality is violated. However, the vi- 
olation never occurs for the product states. This implies that two classical fields in the 
mode-entangled states, although spatially separated, present a nonlocal correlation, which 
means the mode-entangled states can be really regarded as the classical simulation of quan- 
tum entangled states. The relevant experiments are necessary because they can prove the 
feasibility of our scheme, moreover, they might open new perspectives for optical quantum 
computation and communication. 
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Fig. 1: Scheme of the numerical simulation. 

Fig. 2: BPM simulation result for the scheme shown in Fig. 1 corresponding to A = 0. 

Fig. 3: The relation between the intensity differences and the phase modulators: (a) 
A c (0J, A) and 6 U (b) B t (0 2 , A) and 2 . 

Fig. 4: The correlation functions S(9i,9 2 ) for mode-entangled states: (a) \&f), (b) |<£i~)> 
(c) |*+), and (d) |*r). 

Fig. 5: The maximum of \B\ versus the refractive index of the phase modulator 3 . 
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Fig. 1. Scheme of the numerical simulation. 



Fig. 2. BPM simulation result for the scheme shown in Fig. 1 corresponding to A = . 




Fig. 3. The relation between the intensity differences and the phase modulators, (a) 

A c (e x ,X) and 0, , (b) B t {6 2 ,X) and 6 1 . 
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Fig. 4. The correlation functions S{6 1 ,6 2 ) mode-entangled states: (a) |*i }, (b) |<J> t ^ 
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Fig. 5. \b\ versus the refractive index of the phase modulator 3 



